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LINEARIZED COMPRESSIBLE-FLOW THEORY 

FOR SONIC FLIGHT SPEEDS 

By Max. A. Heaslet, Harvard Lomax, and 
John R. Spreiter 


SUMMARY 

The partial differential equation for the perturbation velocity 
potential is examined for free— stream Mach numbers close to and equal 
to one. It is found that, under the assumptions of linearized theory, 
solutions can be found consistent with the theory for lifting— surface 
problems both in stationary three-dimensional flow and in unsteady 
two-dimensional flow. Several examples are solved including a three- 
dimensional swept— back wing and a two-dimensional harmonically 
oscillating wing, both for a free— stream Mach number equal to one. 


INTRODUCTION 

Much of the recent progress in the theoretical analysis of 
compressible— flow fields is attributable to the successful applica- 
tion of linearization methods. Although the basic assumptions used 
in conventional linearized theory appear at first glance to be highly 
restrictive, it has been found that, just as in the analogous case of 
thin-airfoil theory for incompressible flow, the methods have many 
fields of utilization adequate for most engineering purposes. Since 
the basic methods are so well known and depend on such relatively 
simple mathematical tools, it appears obvious that the range of 
applicability of the theory should be explored completely. Such is 
the purpose of the present report. It has been more or less tacitly 
presumed in the past that such applications cannot treat cases for 
which the flight velocity is near the speed of sound. In the study 
of two-dimensional steady— state problems in airfoil theory, this 
presumption is certainly true. The Prandtl-Glauert and Ackeret rules 
for variation of pressure coefficient with free— stream Mach number 
in the subsonic and supersonic regimes, respectively, are clearly 
invalid for Mach numbers near one, since perturbation velocities 
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become arbitrarily large. In this case, linearized theory therefore 
predicts its inability to treat such problems. On the other hand, 
if linear methods are applied to nonstationary two-dimensional 
airfoil and particular steady— state, three-dimensional, lifting- 
surface problems at sonic speeds, a consistent theory results since 
solutions are found which yield perturbation velocities of the same 
order of magnitude as those calculated for free— stream Mach numbers 
of, say, 0o6 or 1.5 » 

Unfortunately, arbitrary thickness distributions at sonic speeds 
cannot be studied by linear theory in the steady state since, in 
general, the theory predicts infinite pressure differences between 
the wing surface and infinity. In the particular case of a yawed, 
symmetrical wing of infinite aspect ratio, the results are, however, 
again consistent with the theory and yield pressure distributions 
which are the same as those determined by using only the component 
of free— stream velocity normal to the leading edge. The derivation 
of this latter result for a free— stream Mach number of one will be 
given. 

The difficulty of not being able to include thickness effects 
in general, together with the uncertainty of the magnitude of the 
viscous effects, leaves the question as to the limitations of such 
a linear theory in application to practical wing shapes. Such a 
question can certainly not be resolved by mathematical reasoning 
alone. The extent to which the fluid medium can be idealized at 
these speeds is left, for the time being, unsettled and it remains 
for experiment to determine whether the consistent mathematical 
results which are obtained from the linearized equations provide 
reasonably exact predictions. In this connection, it should be 
mentioned that the few experimental results available for the total 
lift on thin triangular wings at Mach numbers near one tend to 
confirm the theory. But even if more detailed experimental results 
indicate that further refinements are necessary, there is still 
little doubt but that the linear potential solutions will provide a 
valuable basis for more exact extensions of theory. 

The present report is divided into three parts. In the first 
part, the linearization of the partial differential equation for 
the velocity potential is carried out in some detail for steady- 
state conditions. A by-product of this derivation is the nonlinear 
form of the equation for two-dimensional flow which was used by 
von Karman (reference l) to determine his similarity rules for tran- 
sonic flow. The equation for unsteady two-dimensional flow based on 
the same assumptions is also given. The second part of the report 
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is restricted to two-dimensional unsteady problems for values of 
Mach number near one. The principal contribution of this section is 
the evaluation of the change with time of the pressure distribution 
over an airfoil starting suddenly from rest at a speed close to that 
of sound. Such an idealized problem involves a step function in 
velocity in which the airfoil has zero velocity for all negative 
values and near sonic velocity for all positive values of time. From 
these results the initial build-up of lift can be calculated for Mach 
numbers near one, although the eventual value of the lift cannot be 
found by linear methods. Further application can also be made to 
problems in flutter and gust loads. The third part of the report 
treats the steady— state three-dimensional problem. Both lifting 
surfaces and symmetrical nonlifting wings are considered and it is 
seen that in the former case consistent solutions are obtained by 
particularly simple means , These solutions represent the limiting 
case of both subsonic and supersonic lifting— surface theory and 
give, for example, the same value of lift-curve slope at the speed 
of sound that was obtained for the supersonic triangular wing by 
Stewart (reference 2). 

A list of symbols is given in the appendix. 


PART I - THE LINEARIZED EQUATIONS OF MOTION 
Steady State 

The nonlinear partial differential equation satisfied by the 
velocity potential <j> of an isentropic flow field can be expressed 
in the form 


$ yz »y$z — ®zx — — $xy ^x^y = 0 


( 1 ) 


where the subscript notation is used to indicate differentiation and 
a is the local speed of sound given by the relation 



( 7 - 1 ) 

2 


Mq 2 



( 2 ) 


In this latter equation V 0 and Mq are, respectively, velocity 
and Mach number of the free stream, 7 is the ratio of specific 


4 


NACA TN No. l8?4 


heats (for air, 7=1.4), and V is local velocity. 

Introducing the perturbation velocity potential cp, where 

cp = -V 0 x + $ (3) 

it is possible to express equation (l) in terms of the derivatives 
of cp and the parameters M Q and V Q . To begin the linearization 
of the resulting equation, the coefficients of the second ordered 
derivatives of $ are expanded in Maclaurin series with ascending 

powers of — , — , — . The convergence is assured provided 
Vo Vo 


(7-1) 

2 



+ 


u 2 +v 2 +w 2 ^ 

/ 

v o 


<1 


w 


or, in a slightly modified form, provided 


i , Oo 2 

V 2 -V Q 2 <tfo_ = 5a 0 2 (5) 

O I y_l u 

If the assumption is now made that r—, « 1 so that 

''o *o *o 

second and higher powers in the perturbation velocities can be 
neglected in comparison with one, the partial differential equation 
can be simplified to the form 


c Pxx 1 - M 0 2 


1 + p- + (7-1) M q 2 -p 1 |+ <Pyy + 'Pzz 


- 2 <Pxz -H- Mq 2 - 2 Txy -2- Mq 2 = 0 


(6) 


From this equation all the succeeding expressions will be derived. 

Two— and t hree-dimensional linear equations, Mq ^ 1. — Since 

equation (6) is obviously nonlinear, additional assumptions must be 
made to reduce it to a linear form. Clearly, these assumptions must 
involve the relative magnitudes of all the terms in order to determine 
which ones may be neglected. Perhaps one of the least restrictive 
set of conditions is that: 
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(a) The ratios of the perturbation velocities to 
the free— stream velocity are small enough to be 
neglected when compared to one. 

(b) The velocity gradients at a given point of 
the flow field are all of similar magnitude. 

With the aid of these assumptions, it follows that, to the 
order of the approximations made, the perturbation velocity potential 
Cp satifies the well-known linear equation 


(1 - M 0 2 ) <? xx + <Pyy + <P ZZ = 0 (7) 

In the case of two-dimensional flow, the equation is independent of 
y and thus may be written in the form 

(1 - Mo 2 ) 9^ + q> z2 = 0 ( 8 ) 


Two- and three-dimensional nonlinear equations, Mr> = 1. - 

The study of equation (8) in both subsonic and supersonic flow has 
shown that ior arbitrary lifting surfaces or symmetrical nonlifting 
airfoils the value of the induced velocity u on the surface of a 

fixed geometric configuration is proportional to ([ 1 — M,-, 2 | )”" 1/ ^. 

In all airfoil problems, the value of u becomes infinitely large 
as Mq approaches one, either from above or below, and the basic 
assumptions are thus violated. Such a difficulty led Oswatitsch 
and Wieghardt (reference 3) and Sauer (reference 4) to abandon the 
restriction of linearity and to seek a more exact equation at M = 1. 
Retaining the assumptions underlying equation (6) and setting 
V Q = a* where a* is the critical speed of sound, it follows that 
at M q = 1 the perturbation velocity potential satisfies the equation 


(7+1) 

a* 


<Px 


'Pxx C P zz + ~Z 


xz 


= 0 


(9) 


Since cp x is much larger than cp z as the Mach number approaches 
one, equation (9) may be further simplified to 


(7+1) 


cp cp - cp = 0 

X ^xx ^zz 


( 10 ) 


If, in three dimensions, the perturbation velocities do not 
remain small, equation (6) again supplies the necessary form of the 
differential equation at M c = 1. From the relation V Q = a*, the 
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required expression is 


(7+1) 

a* 


CD CD — CD — CD 

v x Y xx Y yy Y zz 




*7 V = 0 


(id 


Two— and three-dimensional linear equations, M Q = 1* — 

Equation (10) has been used by von Karman (reference l) to establish 
similarity rules for two-dimensional transonic flow and is the basis 
for work continuing at the present time. (See also reference 5») 

If at M Q = 1 the assumptions made in the linearization process still 
hold, it follows from either equation (10) or (8) that the differential 
equation reduces to the form 


<Pzz = 0 (12) 

It is possible, however, to predict independently from this relation 
that linearized methods cannot be applied to the calculation of 
arbitrary airfoil pressure distributions. The range of applicability 
of such an equation is thus almost nonexistent. On the other hand, 
the linearized form of equation (ll) or ( 7) at M Q = 1 is 

V + V 50 (13) 

and from this equation a class of nontrivial solutions can he 
obtained for particular boundary conditions. Both equations are 
of parabolic form in the number of dimensions for which they are 
defined. In the present report, formal solutions satisfying the 
imposed conditions will be obtained in three dimensions for flat 
lifting surfaces with swepi>-back leading edges and for an infinitely 
long, symmetrical, sweptr-hack wing. 


Unsteady State 


The derivation of the steady— state equations for the velocity 
potential was developed in some detail because of the various results 
to be obtained. Similar methods can be used when unsteady conditions 
are to be considered, the differential equation for the velocity 
potential being now in the form 


-L 

a 2 

( ®t*t» + 

+ ( 

>-s; 

2 


a 2 

*yz *y *z 


+ ® 

2_ 

' a 2 


yy 


r zx 


+ 2®y Gytf + 2$) z 
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where t* represents time. The details of the derivation can, 
however, be avoided by referring directly to the equation satisfied 
by the velocity potential for the propagation of sound waves of 
small amplitude. (See reference 6, p. 492.) In this form of the 
equation the Cartesian coordinate system x, y, z is assumed fixed 
in the medium so that free— stream velocity is zero, while the wing, 
which moves in the direction of the negative x axis with velocity 
V Q , generates small pressure disturbances. As a consequence, the 
velocity potential of the field satisfies the well-known wave equation 
in three space dimensions: 

h ‘Pt't* - <*>xx - <Pyy - <f>zz = 0 (15) 

Equation (15) is reducible to canonical form by means of the 
relation 


t — a Q t ^ 

and the three-dimensional form of the equation is therefore 

<Ptt " Txx -< Pyy ~^zz = 0 ( l6 ) 

while in the two-dimensional case independence with respect to y 
yields 

*tt - V - \z * 0 < 17) 


PAST II. TWO-DIMENSIONAL LINEAB PROBLEMS FOE Mo NEAE ODE 
Unsteady State , Mq = 1 

It vas pointed out in the derivation of equation (12) that the 
linear equation for the velocity potential is not applicable to air- 
foil problems in either the subsonic or supersonic regimes for Mq 
near one. The possibility still remains, however, of analyzing 
unsteady flows during the period in which the perturbation velocities 
remain small. As an example of such a problem, consider the case of 
a flat lifting surface at a small angle of attack a starting from 
rest at a velocity Vq near the speed of sound. The perturbation 
potential for such a motion is equivalent to the change in potential 
brought about by an abrupt change a in angle of attack of an air- 
foil flying in a steady— state condition at velocity equal to V Q . 
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The change in load distribution brought about by this maneuver is 
the so-called indie ial load distribution and a knowledge of such 
indicial functions is important for applications using operational 
calculus . 


Figures l(a) and 1(b) 
furnish an insight into 
the nature of the boundary 
conditions for airfoils 
traveling, respectively, 
at supersonic and subsonic 
speeds. The chord is ini- 
tially on the x axis with 
leading edge at the origin 
and trailing edge at 
x = c 0 . With increasing 
time, the wing section 
travels in the negative 
x direction and sweeps 
out a portion of the xt 
plane as shown in the 
figures ( indicated by 
shaded areas). Through- 
out this part of the 
plane the boundary condi- 
tions require that the 
induced vertical velocity 
is — V 0 cx, while elsewhere 
on the plane the induced 
velocities are continuous 
functions of z. In fig- 
ure 2, sketches of the 

airfoil in the supersonic and subsonic cases are shown together 
with indications of the manner in which the disturbance field 
spreads. These wings are presented in xyz space with time as 
a parameter so that their coordinate system is not to be confused 
with that of figure 1. The airfoils are traveling from right to 
left at Mach numbers of 0.8 in the subsonic case and 1.2 in the 
supersonic case, and for a time corresponding to that required 
for the wing to travel a distance of one— third chord length. At 
t = 0 cylindrical waves are induced at each disturbance point, that 
is, at each point of the chord. These waves expand radially at the 
center of the expanding waves moves relative to the initial disturb- 
ance point. At a given instant in time the entire disturbance region 


characteristic traces 



(a) Supersonic wing. 



(b) Subsonic wing. 

Figure 1.— Boundary conditions for 
two-dimensional unsteady— lift 
problem. 


NACA TN No. 1824 


9 


of the airfoil is contained within the closed surfaces shown in the 
figures, the outer surfaces corresponding to the largest values of 


/ 



/ 


Figure 2o- Sketch showing extent of disturbance fields after 
travel of one— third chord length. 

time. In the supersonic case, the pressure distribution over the 
wing reaches a steady— state value as soon as the wing moves ahead of 
the expanding cylindrical wave produced at t = 0 by the leading 
edgeo In the subsonic case, the wing never leaves the disturbance 
field of the cylinders and, as will be seen later, the steady— state 
pressure distribution is approached asymptotically. 

It is apparent from equation (17) that the characteristic cones 
have semivertex angles equal to 45 and that the cones with vertices 
on the xt plane have traces with slopes equal to ±1. These cones 
determine the upstream boundary of the field of influence of the 
vertex point and their cross sections in the plane t = constant 
are the disturbance regions of the cylindrical waves arising at the 
vertex. Thus, perturbations in pressure produced initially at the 
leading edge of the wing section are confined at later time, in the 
modified coordinate system, to the cone with vertex at the origin 
and traces x = ±t. 
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The solution of equation (17) for boundary values of the type 
under discussion has been indicated in reference 7 through considera- 
tion of an analogue problem in supersonic lifting— surface theory. 
Thus, the shaded areas in figures 1(a) and 1(b) are thought of as 
swept— forward lifting surfaces situated in a stream directed along 
the positive t axis at a Mach number Mq = The boundary values 

remain the same; that is, <?z = w = — V G a on the wing and <Pt> <Py, 9z 

are continuous functions of z elsewhere in the xt plane. In 
lifting— surface terminology, the unsteady case for supersonic speed 
becomes a wing with supersonic leading edge, while the case indicated 
in figure 1(b) involves a subsonic leading edge. 

The solution for the wing traveling at supersonic speed has been 
given in reference 7 in a form valid for all Mach numbers greater 
than or equal to one. The expressions for load coefficient Ap/q, 
where 

AP = Pz~P u 

0. ip y ^ 

2 p o v o 


differ analytically in various regions of the xt plane. These 
expressions are: 

Region A (between lines x = -M 0 t, x = — t) 


Ap 4a 



(18a) 


Region B (between lines x =-t, x = t, and x = c 0 -M Q t) 


AR = 
q 


4a 

V^o 2 -l - 


— arc cos 

n 


MpX+t 

x+M 0 t 


+ 


7'Mq^-I /tt + 


JtM 




arc sin 



(18b) 


Region C (between lines x = t, t = 0, and x = c 0 -Mot) 

Ap _ 4a 

q = M 0 


(18c) 


From the pressure distributions it is possible to calculate the 
indicial lift coefficient Cj^(t) as a function of Mq and t. 
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Since 


v*> - 1 


'O 

the following results are obtained: 
First time interval 0 < t < - -0 - 

1+Mq 


Ap 

qa 


dx 



(19a) 


Second time interval 


_£a_ 

1+Mq 


<t< 


C Q 

Mo-1 


- 7 


i /« 


— . . c Q — M t\ n t+M c n -tM. 2 

— - + arc sin ° ) + - lJL arc cos £-2 &- 

L M oV ‘ ' Jt^-1 c o 




(19b) 


Third time interval — 2_<t 

Mq-1 


c ^ (t) =7^ 


(19c) 


These results have been discussed in reference 7 for values of 
Mq greater than one. They still hold, however, for sonic flight 
speeds and, in fact, can be reduced to the expressions: 

First time interval 0 <t< ^ 

2 
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C^U) = k 

Q 

Second time interval — <t 

2 


(20a) 




(t) = «(i + arc Sin£ t^ 



(20b) 


The indie ial lift coefficient is seen to be constant and equal to 4a 

up to tbe time t* = -2- or up to the time required to travel one-half 

2*o 

chord. Following this first time interval the indicial function rises 
monotonically, reaching an infinitely large value as time increases. 



the steady-state load coefficient becomes infinitely large in linear 
theory for 1^=1. This means that the theory cannot be used to predict 
the complete extent of the C^(t) variation with time but that 

during the earlier part of the motion the assumptions remain valid. 


2MJ 


Co 



Figure 3 •— Indicial— lift-curve 
slope for Mach numbers between 
0 and 1.4 shown to time required 
to travel 12 half— chord lengths. 


In figure curves of Ct (t) 

2M c t 

are plotted as functions of — y — 

c o 

for values of Mq= 1 as given by 
equations ( 20a ) , ( 20b ) , and 
M q =1.2, M q = 1.4 as given by equations 

(19a), (19b), and (19c). Also 
included in the figure are variations 
of Ct (t) for M Q =0 as calculated 

sc, 

from Wagner's results (reference 8) 
by R. T. Jones (reference 9) and- 
also for Mq= 0.8. The derivation 
of results leading to the Mg =0.8 
curve will be given subsequently 
in this paper. The value of Cj^(t) 

at Mq= 0.4 for a short interval 
of time is also drawn. The dashed 
portions of the curves were not 
calculated but were drawn to agree 
with the known asymptotic value 
of the lift function. 


NACA TN No. 1&24 


13 


Application of the indie ial lift function at Mq= 1.- Once the 


indicial lift function is known, it is possible to determine the 

lift corresponding to a given variable motion. Consider, as an 

idealized example, the case 

where the airfoil experiences 

an abrupt rising and sinking 

motion at regular time intei>- a(t') 

vals . Such a motion involves 

abrupt plus and minus angles 1 j j 

of attack without rotation or — — 1 1 \ 1 / 


pitching so that a(t f ) is 
given by the meander or square- 
wave function shown in figure 
4(a). In this example the 
variation of a is such that 
the curve for C L (t f ) can be 

calculated easily. In figure 
4(b) C L (t f ) is shown for the 
case in which the discontin- 
uities occur at intervals of 
time equal to c 0 /V 0 , that is, 
after each chord length of 
travel. The principal point 
of interest in this example 
is the fact that such a 
motion yields no excessive 
value of lift or perturbation 
velocities and the entire 
analysis is within the frame- 
work of linear methods. 


i 

i 


(a) Impressed angle of attack. 

cyo 



(b) Resulting variation of lift. 

Figure 4.— Lift resulting from 
square-wave angle— of -attack 
variation. 


When the variable motion is more complex in character the lift 
coefficient can be expressed by means of Duhamel’s integral. Corre- 
sponding to the angle— of-attack variation a(t’) as a function of 
time, the lift coefficient C^t’) is given by the expression 


c L (t*) 


aF 


Cljt'-T*) a (t») dr* 


( 21 ) 


In analysis related to equation (21) it is convenient to employ 
techniques associated with the use of the Laplace transformation. (See 
reference 10.) Thus, if the Laplace transform f(s) of the function 
f(t) is defined by the relation 
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f ( s ) = / e - ® 1 f(t) dt 


then equation (2l) can "be rewritten in the form 


Cl(s) = sCi^(s)a(s) 


( 22 ) 


Consider now the case of a lifting flat plate oscillating 

harmonically without pitching at a frequency <n and maximum angle 

of attack equal to a . Setting 

max 


/ , \ lout iaua 0 t 1 

a(t) = °4nax ® = “max e ° 


(23) 


then equation (22) yields 


_ b £o 

— 2 - erf /ST 

Sax s-ico \*lc 0 xj^ V 2/ 


Cl(s) = 



(24) 


By straightforward manipulation, the inverse transformation of 
equation (24) can he shown to give 


SL. . ke M -ll ♦ /4 e" lv 


.11 


C(oDit— v) — i S(o)tr- V) 


-2 Ji Ni(ort, v) - N 2 (cot, v)| - | arc (25) 


where 


v = ^o _ tu'cp 
2 " 2T~ 
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a)’ is true impressed frequency 
ioi't 1 

( a = cu max e )• 

C(aot— v), S(cot— V) are Fresnel's 
integrals (reference 11). 

N ]_ ( cot , v) , N 2 (a>t , v) are integrals 
defined in the appendix. 

If the response to a cosine 
variation of a is required, 
only the real part of equation 
(25) is used. Such a response 
in the early stages of the 
maneuver is shown in figure 5. 

For very large values of time 
C L as given hy equation (25) 

approaches the value 



Figure 5.— Lift resulting from 
cosine-wave angle— of-attack 
variation. 


CL — iot f r 1 iv \ 

!C’* e \vI7v e -erf/TvJ (26) 

from which both the amplitude and phase shift of C L resulting from 
either an impressed sine or cosine variation of a can be readily 



Figure 6.— Amplitude of oscillatory lift resulting from a cosine 
angle-of-attack oscillation (without pitching) at M Q = 1. 
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determined. Figure 6 shows the amplitude of lift oscillation corre- 
sponding to continuous angle— of— attack oscillation, as determined 



Figure 7« — Variation of cos 03t 
with t for various values of 

03. 


from equation (26), plotted as 
a function of v = It is 

apparent that as 03 approaches 
zero, approaches infinity; 

that is, as the impressed wave 
approaches the "step" function 
(fig. 7)* the lift coefficient 
approaches infinity. This 
result is in agreement with 
equation (20b) as t approaches 
infinity. As the frequency 
parameter V is increased, 
however, the value of C j,^ / a max 


is reduced and reaches a minimum of about 3*4 for a value of V =0.9* 
For a speed of sound around 1000 feet per second and a wing chord of 
6 feet, this would correspond to a frequency of 47*7 cycles per second, 
a value well within the range of practical flutter frequencies. It is 
interesting to note that figure 6 also shows that as the frequency of 
oscillation becomes large (i.e., v>3) the value of CL^ y /°4mx 

approaches the value 4. This is the same as the value for Cj^(t) 

in the early stages following a step variation of a. 


Unsteady State, M Q < 1 

It has been pointed out that, in the determination of the 
indicial lift function for a wing traveling at subsonic speeds, the 
lifting— surface analogue involves the calculation of load distribution 
over a swept— forward wing with subsonic edges. This means that 
recourse cannot be made in the solution to the simple source distri- 
bution method used in reference 7 to treat the 1 case. Since, 

however, a portion of the leading edge in the present case is still 
supersonic, the problem is particularly adapted to lifting-surface 
methods developed by Eward in reference 12. Figure 8 indicates, as 
in figure l(b), the geometry associated with the boundary conditions. 

In the Eward analysis the solutions, as in the previous case, are 
calculated for various regions. As an aid in identifying the different 
results the sketch denotes these regions by Roman numerals. 
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The load distributions 
and indicial lift functions 
will be given here for values 
of t up to the time when 
the characteristic from the 
trailing edge first crosses 
the leading-edge trace, that 
is, for 

1-Mq 

or 


0 - t* = 


a o ( 1 “^o ) 



Figure 8.- Regions used in the 

study of subsonic unsteady lift. 


This period covers the time 
in which the wing travels 
2M 0 /(l-M Q ) half chords and, 

since the present analysis 

is concerned with values of Mach number near one, will in some cases 
extend beyond the range of the linear theory. 


The following results are obtained for load coefficient: 
Region I (between lines x=t, t=0, anri x=c Q — t) 

Ap 4a 

<1 = Mo 

Region II (between lines x^M^, x=t, and x=c Q — t) 

Ap _ 8a / Mp j~ 
q _ nM 0 U+Mo/ * 


(27a) 


t— x 

M^+x 


+ arc tan 


m 


(27b) 


Region III (between lines x=c Q — t, x=t, x=— t + 

A 


2c, 


1+Mq 


and 


i-^y 


t 
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Ap 

q 


8a ] t— x 

jt ( 1 +Mq ) J x+Mot + 


4a 

kMq 


arc sin 


2x-t(l-^) 

(l+Mo)t 


+ arc sin 


2(c 0 -x)-t(l+M 0 ) 

t(l-Mo) 


(27c) 


Region IV (between lines x=t , x=c 0 -M Q t, and x=c Q — t) 


Ap _ 8a 
q kMq 


arc sm 


/x-c Q +M Q t 

t(M 0 -l) 


(27d) 


Region V ^between lines 


2c n 

x= — t + — — 


1+Mq' 


x=c 0 — M 0 t, and t 


c o 

1-Mq/ 


Ap 


l6a 

« 2 (l+Mo) 



| - ER^k') - KE(i,k») + KF(t,k’) 


2a , x 4a , 2x— t(l— M 0 ) 32a K / c -x-M t 

vT a^c sin — + —rr- arc sin .. - vt — + -5- / - — y 

t kM o (l+M^t « 2 ( 1 +Mq) v (l-^o 2 ) (x+t) 


2a 
* Mo 


arc sin 


2c 0 -t(l+M 0 ) 

t(l+Mo) 


+ N 


3 


where 


k* =./l-k 2 


k 



2c 0 

(t+x) (1+Mo) 


(27e) 



i|r = arc sin 
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V- ZSL 


c n </2 


N 0 = 


ka 


Mo 


-J l (v-Vi)(u+V 1 )- 


X 


| sin -i (u+-7i)(l-^o)[(v-^i)(l+Mo)-c 0< /2 3 +2vi [ v( 1+Mp )— c 0 \/^-u( 1 -Mq ) ] ^ 
^ [v(l+Mo)-c 0 ^ + Vid-^)] [ u( 1 -Mq ) - vi(l+M 0 )] J 


dvi 


where 


u = — (t-x) 


v = -7= (t+x) 

E^k) J* incomplete elliptic integrals 


K, E complete elliptic integrals 


In figure 9(a) the growth of pressure distribution with time is 
shown at subsonic speed for the period of time covered by equations (27), 
For purposes of comparison, 
pressure changes calculated 
from equations (18) are shown 
in figure 9(b) for supersonic 
flight velocities. 

Equations (27) suffice 
to determine the initial 
growth of indie ial lift coeffi- 
cient at subsonic speeds. 

Such results were given in 
figure 3 at M Q =0.8 along 
with the calculated growth 
for about one chord length 

The 

value of C T ( t ) at t = 0 



of travel at M Q =0.4. 
% (t) 

is, for all flight speeds, 
equal to b/V^. 


Figure 9*- Pressure distribution on 
wings receiving sudden angle-of— 
attack change at t = 0. 


Expressions for Cj^(t ) are as follows: 
First time interval 0 <t< - C 0 

1+Mo 
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C Iu (t) 


4 


Cp-fc(l-Mp) 


Second time interval 


1+M, 


<t < 


1-Mo 


(28a) 


Cl^t) 


_4_ r+c 0 -3t(l^) / g c„-t(l-M/) 

*Vo l itM c y 


+ * 



+ 1+3M q 

1+Mq/ ( 1+Mq ) 2 


\^2t ( 1 +Mq )— 2c 0 1 [2c 0 -t(l-^o 2 )] 


+ ^o^td+Mp) 

1+M 0 


arc sin 


t(l+M Q )-c Q 

t(l+Mo) 



+ [2c c — td+Mo)] arc tan 


* 


, c c _M o t 


t(l+Mo)-c 0 J C ° 2c 


&) dx 
qa> 


1+Mp 


-t 


V 


where gi- ven ^y equation (27e). 


(28b) 


PART III - THREE-DIMENSIONAL LINEAR PROBLEMS FOR Mp NEAR ONE 


Steady State 

General solutions for arbitrary Mach numbers .— Two methods of 

attack are available for the solution of linearized problems at sonic 
speeds. In the first place, solutions to equation (13) can be 
written formally and the extent to which these solutions satisfy the 
original assumption can then be investigated. In the second place, 
general solutions of equation (7) can be studied in the limit as M Q 
approaches 1. Since this latter method furnishes added information 
concerning the variation of the variables with Mp, it will be used 
first. 
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In linearized theory the boundary— value problems of wing theory 
are concerned with two separate properties of the wing: the thickness 
effects and the effects produced by the twist, camber, and angle of 
attack. The first is called the nonliving case and the second is 
the lifting case. Solutions of equation ( 7 ) for 1 are given 

in reference ( 13 ) as follows: 

In the nonlifting case 


<P(x,y,z) 



Aw 0 (x 1 ,y 1 )d^: 1 dy 1 
^(x— xi) 2 — p 2 (y-yi) 2 — P 2 z s 


(29) 


where P “V Hd - 1 and Aw 0 = 2w Q where w c is the vertical pertur- 
bation velocity on the wing and therefore related directly to the 
slope of the wing surface relative to the x axis. The integration 
region t is the area on the wing within the Mach forecone from the 
point x, y, z. 


In the lifting case 


T(x,y,z) 


_1 

2it 


p p 3 g z A(p g) (x 1 ,y 1 )dxidy 1 
'l [ ( x-xi ) 2 -p 2 ( y-y 1 ) 2 -P 2 z 2 ] 3/ ^ 2 


(30) 


where A<p Q is the jump in the val ue of t he velocity potential in 
the plane of the wing. The sign | denotes "finite part" of 

the integral and introduces special integration techniques. (See 
reference 13 . ) 

Equation (29) expresses the velocity potential for the symmetri- 
cal wing in terms of an integral involving supersonic source distri- 
butions while equation (30) employs doublet distributions. In the 
two cases the distributions are determined from the geometry and the 
load distribution over the wing, respectively. 

Source and doublet distribution effectiveness at infinity .— It 

is well known that the lift, drag, and pitching moment of a given wing 
may be calculated either from direct integration of the local pressures 
on the wing or by means of momentum considerations where the induced 
velocities of the wing are determined at an infinite distance and 
the desired forces are related to an integration over a control surface 
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enclosing the wing. In the three following sections the latter 
approach will he considered and the limiting value of drag at M Q =1 
computed. The initial portion of this theory requires the evaluation 
of source and doublet effectiveness at infinity and the concept of 
equivalent source position, an idea which appears to have been given 
first by W. D. Hayes in reference 14. 

Consider, as in figure 10, a point P with coordinates x, y, z 
lying within the induced field of a supersonic wing. The Mach f ore- 
cone from P is given by the relation 


x— Xi 


= PV^y-yi ) 1 2 + (z-zi ) 2 


(31) 


where x±, zi are running coordinates of a point on the surface 

of the cone. Introducing polar 
coordinates 



Figure 10.— Coordinates used in 
study of supersonic source. 


y = r cos 6, z = r sin 0 

and rewriting the abscissa of 
P in the form 

x = x 0 +pr 

it follows that the trace of 
the fore cone in the zi =0 
plane is, in the limit as r 
appr oache s infinity , 

Xi = pyi cos 9 + x c ( 32 ) 


It is, moreover, possible to show that the effect on the velocity 
potential at the point P as r approaches infinity is the same 
for all points (xi, y 1? o) for which x^.— Pyicos 0 = constant. The 
value of this effect is 


cp = 


1 

2* v/2(3r(x 0 — xx+Pyxcos 0) 


(33) 


and follows from the asymptotic evaluation of the supersonic source 
potential 


1 

2 n J l x-xi ) 2 -P 2 ( y-yi ) 2 -(3 z 2 


<P 
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for large values of r. The potential at P for the source at 
( x i> 7i> 0) thus the same as for the source shifted along the 
trace to (x 1 -f3y 1 cos 6, 0,0), the intercept of the trace on the x 
axis. For Mach numbers near one, equation (33) can be rewritten 


cp = ^ 

2 *V 2 P r ( x °- x i)" (3*0 

and is equivalent to the potential at P for the source at (x;l, 0, 0). 
The induced velocities at P due to a source at (x ly y lf 0) follow 
immediately, for arbitrary M 0 and for M 0 near one, from the 
gradients of cp in equations (33) and (3^)* It is important to note 
that equation (33) is a function of the azimuthal angle of P so 
that, in general, a source does not have a fixed equivalent position 
with respect to its potential at infinity; equation (3*0> however, 
is independent of the azimuth 0. 

The source— sink potential is applicable to the study of 
symmetrical nonlifting wings. When lifting surfaces are to be analyzed, 
the doublet potential 


2it [ ( x-xi ) 2 -0 2 ( y-yi ) 2 -p 2 z 2 ] 3/2 

% 

must be considered and the question of equivalent doublet position 
with respect to the potential at infinity arises. In this case the 
doublet position can again be shifted parallel to the trace of the 
Mach cone from P at infinity and the potential at P is given by 
the expression 


cp = 


P 2 z 


2it[2|Br(xo-^xi+PyiCOs e)] 3 ^' 
and, for Mach numbers near one. 


(35) 


cp = 


3 2 z 


2jc [ 2Pr ( x 0 — X! ) ] 3 / 2 


(36) 


Momentum relations .— The vectorial force ^ on an aerodynamic 
body inside a control surface S is given by the surface integral 
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dS 






where vector notation is used and 


(37) 


o subscript indicating free— stream condition 

p,p local static pressure and density 

V local perturbation velocity vector 

For the purposes of the present report, equation (37) will be 
modified according to the assumptions of linearized theory and the 
surface S restricted to a semi— infinite circular cylinder of 
radius r, its axis of symmetry lying along the x axis, and with one 
face in the x=0 plane while the other face is at x=constant. (See 
figure 11. ) 


From linearized theory, 


_P_ 

Po 


1 - v 


u 

Vo 


and 


P - Po = 


V q u + |(u 2 +v 2 +w 2 ) 


+ 5 »o M o 


2 U 2 



The end faces of the cylinder 
may he denoted, as in the 
figure, by I, II, and the 
curved surface hy III. Then in 
supersonic flow, if a distribur- 
tion of sources is restricted 
to a region downstream of I, 
the drag D on the hody corre- 
sponding to the source distribu- 
tion is given hy the expression 
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p o 


f J (M 0 2 -1)u 2 +v 2 +w 2 J 


dydz— p Q 



uv r rd©cLx 

(38) 


where v r is the radial component of the perturbation velocity. No 
loss in generality results, moreover, if the surface II is moved 
infinitely distant downstream and the radius of the cylinder is made 
arbitrarily large. The notation II and III will henceforth refer to 
this particular configuration. 


If the drag of a lifting surface is to be calculated, the 
surface and its vorticity wake are replaced by doublet distributions 
and in that case the integral over region II in equation (38) is 
called the vortex drag of the body while region III yields the 
wave drag . It has also been shown (see, for example, reference 17) 
that the vortex drag of a supersonic wing is a function only of its 
span load distribution and is equal to the induced drag at subsonic 
speeds for the same span loading. If a finite nonlifting body is 
considered, each of the velocity components in region II is attenuated 
in such a manner that its contribution to the vortex drag is zero. 

The integration over region III again provides the wave drag for the 
nonlifting body. 


The combination of the results given in this and the last 
section provides a method for finding the wave drag of an arbitrary 
body. The first step is the determination of the source— sink or 
doublet distribution corresponding to the body and then, by means of 
the principle of equivalent positions, the sources or doublets are 
moved to the x axis. The wave drag is then calculated from equa- 
tion (38) once the induced velocities on the control surface are 
known. In the next section the wave drag will be written in a 
different form and the drag at sonic speeds will be investigated. 

y This analysis will also provide 
some insight Into the range of 
validity of the sonic theory. 

Evaluation of wave drag as 
Mr. approaches one .— In order to 
study the drag of a symmetrical 
body at zero angle of attack, it 
is convenient to consider the 
general expression for the velocity 
potential given in equation (29). 
Introducing first the transforma- 
tion (fig. 12) 

Figure 12.- System of axes in 
transformation equation (39). 
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where 


t 


= — 7i tan n + xi 


T) = 7i sec [i 


(39) 


tan \i = 3 cos G 


equation (29) becomes 


<pU,y,z) 


i_ r r cos n aw 0 (i ,n)dgdii 
2 * J Jt ( X-Xi ) 2 -P 2 '( y-yi ) 2 -p 2 z 2 


Since, however, it has been shown that cp evaluated infinitely fax 
away from the wing does not change if a source is moved along the 
line |=constant, it follows that the source strengths can be 
integrated along these lines. The second integration is then along 
t]=0 where, from equation ( 39 ) > !=xi and. the value of the potential 
at an infinite distance is 


cpU,y,z) 



cos \i dxi 

V (x-xi) 2 — (3 2 y^— p 2 z 2 


J Av 0 (xi,n)&n 


Setting 


fUi,d) 


it follows that 


— COS |! 



Aw 0 (x!,Ti)dTi 


cp (x,r,0) 



f(xi,d)dxx 
V (x— xx) 2 — 0 2 r 2 


(40) 


and this is the same as the potential for a body of revolution with 
source strength per unit length given by f(xx). The induced velocities 
corresponding to the potential in equation (40) are found to be, after 
first integrating by parts and using the notation d/dxx f(xx,d) = 
f’(xx,|i) together with the relation f(0,ti) = 0, 
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and 




f *(xi,l-i)dxi 
J ( x-Ki ) 2 — (3 2 r 2 


vr 


Sep 

Sr 


-1 

2itr 


X-(3r (x-xi)f t (x 1> |a)dxi 
7(x-x i) 2 -p 2 r 2 


(41) 


(42) 


Asymptotic values of the velocity components for large ‘values of r 
are readily seen to he, after first setting x=x 0 +3r, 


and 


u 



f t ( xi,p)dx i 

/x 0 -xi 



f , (x 1 ,p)dxi 
'Txo-xi 


(43) 


(44) 


Equations (43) and (44) may he used together with equation (38) 
to give for the value of drag the expression 





f t (x 1 , l i)dx 1 p° f»(x 2 ,u)dx 2 

■v/Xq-Xi J q 


(45) 


Assuming that the body is of finite length so that f(x)=0 for 
x>Z reversal of the order of integration yields the relation 



f *( x i^d)f ’(x 2 ,n)2n |xi— x 2 |dx x dx 2 


(46) 
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If equation (46) had "been derived for a body of revolution, then 
f(x) would have been independent of the angle p and in that case 
the expression for drag would reduce to the form 


D = 



f '(xi)f '(x 2 ) 2 n|xi— X 2 |dxidx 2 


(46a) 


/ / 

This expression was given by von Karman in reference 15 . 

For the study of the drag of a lifting surface, consider now the 
general expression for the velocity potential given by equation ( 30 ). 
The doublet distribution occupies in this case both the wing plan form 
and the wake since the jump in 9 exists also in the vortex wake. 

By use of the transformations in equation (39 ) > equation (30) becomes 


9 (x,y,z) 


P 2 z r r cos p z^pU,T))djdT] 

2* Jr J [ ( x— xi ) 2 — 3 2 ( y— yx ) 2 — P 2 z 2 l 3 / 2 


and, exactly as in the case of the source distribution, this can be 
reduced to 


m , v 3 2 z 
9 (x,y,z) = — 


cos p dxi 

[ ( x-xi ) 2 -f3 2 y 2 — 3 2 z 2 ] 3/2 k 


A9(xx,Ti)dTi 


Setting 


6 (^i,p) = cos p 


A9(xi,Ti)dTi 


it follows that 


9 (x,r, 0 ) = 


P 2 r sin 6 
2 * 


l 


X g(xx,p)dX! 

[( x _ Xl )2_p2 r 2]3/i 
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Integrating by parts and using the fact that g(0,p)=0 


and 


u = 


Sep 

Sx 


sin 9 
2«r 


l 


pX—&T p 2 r 2 g t( Xl ^) dXl 


[ ( x— xi ) 2 — 0 2 r 2 ] 3 / E 


(*7) 


_ Sep _ sin 9 
Vv "ot 2?tr 



(x-x 1 )g , (xi,^)dx 1 

y(x- Xl ) 2 -P 2 r 2 


X_f3r (x-xi)3 2 rg t (xi,p)dxi 
[ ( x— xi ) 2 H3 2 r 2 ] 3 / 2 


m 


where g’(xx,|i) indicates g(xx,M-). Setting x = x +0r and 

OXx 

letting r approach infinity, the asymptotic expressions for 
equations (47) and (48) become, if g'(0,p) = 0, 


and 


u 


— sin 6 
2rt 



g* t (xx,p)dx 1 

«/Xq-Xi 


sin 9 „ /T r X ° g'’(xx^)dxi 

‘ 2 ” ■JZ'i 7 ^ 


(49) 


(50) 


The relations just derived may be used in conjunction with 

equation (38) to give the wave drag of a lifting surface. This 
result takes the form 


D = 





g 1 , (x 1 ,q)dx 1 r ° g* > (x g ,p)dxg 
tj X Q -Xi .1 q fj x o X 2 


( 51 ) 
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In the wake of a lifting wing the function g”(xi,|i) = 0 and if, 
moreover, 


J g’ '(xi,n)dxi = g’(c 0 ,h) - g'(0,n) = 0 

o 

reversal of integration in equation (5l) yields the simpler expression 


D = 





g* , (xi,n)g* *(x 2 ,p)ln jxi-x 2 | dxx dx 2 

(52) 


It is possible to draw some general conclusions from equations 
(46) and ( 52 ) regarding the wave drag of wings and bodies of 
revolution without the necessity of detailed applications to particular 
configurations. It is apparent immediately from equation (46a) that 
the wave drag of a body of revolution at zero angle of attack is 
independent of Mach number. This conclusion does not apply, however, 
to the nonlifting wing since the distribution function f(x,p) in 
equation (46) contains the variable p which, in turn, is a function 
of both 0 and (3. As Mo approaches one, the study of the non- 
lifting wing is divided most conveniently into two parts, depending on 
the behavior of f(x,p). 

Consider first the more general situation in which f(x,p) is 
not zero; that is, the case in which the number of sources does not 
equal the number of sinks along the line | = constant. This means, 

when Mq is 1, that an unequal number of sources and sinks appear in 
the transverse or yz plane and, if equation (46) is applied, either 
a finite or an infinite value of drag can result. The limiting value 
of drag at sonic speed, obtained from integrations of surface pressures, 
was given by Stewart and Puckett in reference (16) for several wing 
plan forms, all of which had nonvanishing values of f(x,p). If the 
pressure distribution is calculated, however, the local pressure 
coefficients are seen to become infinitely large as sonic speed is 
reached, even for the body of revolution, so that the assumptions 
of the linear theory are violated and the reliability of the drag 
predicted by equation (46) can in no case be assessed even though 
the predicted values remain finite. Equation (43) shows also that 
when control— surface methods are used to compute drag at Mq = 1, 
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the x component of induced velocity increases indefinitely when 
f(x,p) is not zero and that the theory is, therefore, no longer 
consistent. 


In the very special second case, that is, when f(x,p) vanishes 
for all values of 0, the analysis just presented breaks down at 
equation (40). It is clear, however, that in this case there are 
equal numbers of sources and sinks in the | = constant plane nnH 
the behavior of the flow field at infinity is, therefore, exactly 
the same as that which would have been produced by a distribution 
of doublets. Equations (49) and (50) give the velocities induced 
at infinity by an arbitrary doublet distribution. These induced 
velocity components are, in terms of p, one degree higher than the 
similar components for the nonlifting case. The values of both 
u and v r can thus be expected to approach zero for all values of 
Mo as r approaches infinity for any flow field generated entirely 
by doublets or by an equal number of sources and sinks. It follows 
then that the linearized theory for lifting surfaces (generated 
entirely by doublets) and for bodies with thickness distributions 
such that f(x,|i) vanishes (generated by an equal number of sources 
and sinks in all £ = constant planes) is entirely consistent as 
M q approaches one and, in particular, for Mq equal to one. This 
being true, it follows immediately from equation (52) that the wave 
drag of a lifting system is zero at sonic speed. 

Thickness solutions at Mp = 1 .— A swept-back wing of constant 

chord and infinite aspect ratio is an example of a practical aero- 
dynamic shape for which an equal number of sources and sinks occur 
in every yz plane. (See fig. 13.) Consider the case in which the 
wing cross section is diamond shaped with a slope equal to X in a 
plane normal to the leading edge. Then, in a transverse plane, 
(section BB of fig. 13) w Q equals ±VoX cos >|r, the minus and plus 
signs applying, respectively, to the left and right of the ridge line. 
Accordingly, the solution of the problem can be written in terms 
of a distribution of sources, thus 


-,x cot \|r 


cp = — 

2it 


x — 


2 cos i|r. 


cot 'k 


2m 


x + 


2 cos i|r 


cot 


V Q \ cos i Zn [(y-yi) 2 +z 2 ] dy x 


i|r 

V Q X cos t Zn f(y-yi) 2 +z 2 ] dyx 


( 53 ) 


x cot \(r 
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The value of d<P/dx can immediately be found to be 



from vhich it is apparent that as r = */y 2 +z 2 becomes infinitely 
large, c5<p/c)x approaches zero. In the plane of the airfoil, that 
is, for z = 0, Scp/cht becomes 
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(to) 

I - T o X COS ♦ ; n 

(y tan x— ^c 0 sec V)(y tan x+ ^c 0 sec t) * 


0 « tan 

L (y tan \|r— x) 2 J 


and, using the definition for pressure coefficient, C p = - — , 
this becomes V o 


°P =~ 


2 X cos 'jr 
jr tan \|f 


Zn 




sin t — x cos i|r 



(55) 


Equation (55) can be derived by entirely different methods. 
Perhaps the most direct of these alternative derivations is the one 
introduced by E. T. Jones in reference 18. The general statement used 
in that report is that the component of translational velocity of a 
cylindrical body in the direction of its long axis has no effect on 
the motion of a frictionless fluid. Hence, the pressures over the 
wing shown in figure 13 are the same as those over a wing moving 
normal to a free stream with a velocity V Q cos i|r. Using the Prandtl— 
Glauert correction to the thin airfoil solution of a two-dimensional, 
diamond— shaped, nonlifting section exposed to a free stream with 
velocity V D cos i|r, one obtains, for Mq cos t<l. 


9 = 


r c 0 / 2 

1_ v o 

C0/2 ^l-M 0 2 cos 2 \(r 


Zn[(x'— xi’ ) 2 +z 2 ]dxi* 


(56) 


where w Q is the vertical induced velocity on the upper side of the 
z = 0 plane and x’ is measured normal to the leading edge. If 
this solution is referred to the axial system of figure 13 by the 
transformat ion 


x* = x cos i|r — y sin 

and the integration is performed after taking the partial derivative 
with respect to x, the resultant expression for pressure coefficient 
is 

( l c ° 

Vxcos i|f — ysin \|r 


C p 


2\cos2 \jr 

« yi-Mo 2 cos 2 i|f 



(57) 
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At sonic speed this equation reduces immediately to 


(: 


1 




x cos \jr — y sin 


-1 


(58) 


which is identical to equation (55)* The result expressed by the two 
equations is, of course, not new. The significant point is that the 
same variation in pressure coefficient was obtained by two widely 
different avenues of approach and that the result obtained from the 
particular methods applicable to sonic speed theory is in agreement 
with that derived from more conventional analysis. 

Lifting— surface solutions at Mq = 1 .— It should be mentioned at 

this point that Robinson and Young (reference 19) have shown by means 
of linearized theory that supersonic triangular wings and subsonic 
elliptical wings of the same aspect ratio have values of lift— curve 
slope which approach a common and finite limit as Mq = 1. The 
present section of this report is concerned only with the study of 
lifting surfaces at a fixed sonic velocity but the results to be 
obtained are in agreement with the limiting values of reference 19. 

A further application of the results in this section can be 
made to the case of very low aspect ratio wings at arbitrary Mach 
numbers. This viewpoint of the theory was first presented by 
R. T. Jones in reference 20 and applied to triangular wings while 
in reference 21 extension was made to include pointed wings on 
slender bodies of revolution. This duality of interpretation, that 
is, to all aspect ratios at sonic speed or low aspect ratios at all 
Mach numbers, applies to all solutions of three-dimensional problems 
obtained from equation (13)* In the subsequent analysis, attention 
will be confined to swept— back plan forms of lifting surfaces with 
pointed vertices and thus doublets will be used exclusively. 

In application, the two types of boundary conditions to be 
considered are as follows: 

1. Boundary-value problem of the first kind, loading specif ied .- 

It is given that Au Q = u^— ■ u^ = 0 over the xy plane except for 
the region occupied by the wing where 2x1^ = — 2u^ = Au Q = f(x,y), 
the function being determined by the specified loading. Over all 
of the xy plane, the imposed conditions are Aw Q = 0. 
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2 . Boundary— value problem of the second kind, surface specified. 


Over the xy plane, the imposed conditions are Aw Q = 0 
where and, except for the region occupied by the wing, Auo = 
Over the region occupied by the wing w D = w u = w i = 
f(x,y) is determined by known camber, twist, 

(The delta notation again indicates the jump 


variable at the z = 0 plane . 
Subscripts u and Z indi- 
cate conditions on the upper 
and lower surface, respec- 
tively, of this plane.) 

The nature of the differential 
equation shows that the value of 9 
is a consequence of boundary condi- 
tions along lateral strips. If, as 
in figure lk, the two leading 
edges are given by the expressions 
y = bi(x) and y = bs(x), the 
velocity potential is express- 
ible in the form 


every— 

0 . 

f (x,y) where 
wrirl angle of incidence, 
in the value of the 


y 



Figure l 1 !-.— Swept-back plan 
form with curved trailing 
edge. 


cp( x >y> z ) 


z 

2 * 


r h 2 Aq> n (x,yi)dyi 

ba. (y-yi) 2+z2 


(59) 


If the boundary— value problem is one of the first kind, the 
general expression for cp follows from a direct integration after 
noting that 

f y 

ACp Q (x,y) = I Au Q (x,yi)dyi (60) 


Since, moreover, load coefficient Ap/q is related to Au Q by means 
of the equation 



it follows that the velocity potential 9 can be found for any 
prescribed load distribution of a given plan form. The value of 
vertical induced velocity, evaluated at z = 0, then suffices to 
calculate the twist and angle of attack of the wing. 
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If the boundary— value problem is one of the second kind, the 
vertical induced velocity is given on the wing and the load distri- 
bution is to be found. In this case the use of equation (59) leads 
to the consideration of an integral equation. Since, however, this 
integral equation is a common one in aerodynamic theory, certain 
established methods may be applied to it. 

After noting that A9 Q (x,y) =0 at the leading edge, integra- 
tion by parts and introduction of the relation 


Av, 


dA<P 0 

6y 


yields for perturbation potential the expression 

>b? 


p Up y. y. 

9 = ^ J Av 0 (x,y!)arc tan —— dy x 


(61) 


In the limit as z approaches zero the derivative of 9 with respect 
to z reduces to the form 


_1 r b Av Q (x,y x )dyi 

2* j _ b y-yi 


(62) 


For a given distribution of w D over the plan form of the wing, 
equation (62) represents an integral equation to be solved for 
Av Q (x,y) subject only to the condition that the Kutta— Joukowski 
condition is satisfied at all subsonic trailing edges. Once Av Q (x,y) 
is determined it follows that 


and 


A9 0 = / Av (x,y x )dyi 
— b 

A£ = _2_ 6A9 n 

q V c 8x 


(63) 

(64) 


In the present report the solution to the wing plan form shown 
in figure 15 will be presented. The value of A9 Q which satisfies 
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equation (62) is, in region 1 


A<P = — 2w Jx 2 tan 2 0-y 2 


(65) 


and in region 2 

A<p = — 2 w q x tan 0 


E(V 0 ,k 0 ) — k Q ' 2 F(\|r 0 ,k 0 ) 


where E and F are defined in the appendix and where 


( 66 ) 


,,, . / x 2 tan 0— y 2 

t = arc sm/ — 

,/ x 2 tan 2 0-ax 2 


(67) 


V 


ai 


x tan 0 


= Vi-V 


( 68 ) 


The equation y = a x (x) of 
the trailing edge for which 
equation (66) is valid is 
given by the formula 


a i = 


ko' 


■®o ko 1 Kq 


(69) 


which expresses ai expli- 
citly as a function of 


ai 



This partic— 


Figure 15.— Dimensions and regions 
used in discussion of swept-back 
wings. 


/ 

^x tan 0, 

ular choice of trailing— 
edge shape was used to 
simplify the analysis. The 
resulting plan form approaches a constant-chord wing as the span 
increases. The variation of ai with x is given in figure 16; and 
figure 17 shows the relation between aspect ratio and span. 

The loading coefficient is given in the two regions (defined in 
fig. 15) as follows: 
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c 0 tan# 




1.0 1.4 1.8 2.2 2.6 

s 9 

Co tan# 


Figure 16.— Graph showing trailing— 
edge position of the swept-back 
wings studied. 

Region 1 


Figure 17.— Relation between 
aspect ratio and wing semi- 
span. 


Ap 

qa 


4x tan 2 9 
J x 2 tan 2 0— y 2 


(70) 


Region 2 


= 

qa 


4 tan 0 


E(* 0 > k o> 


x tan 0 


— zfrsif !l° f (♦ ,k ) 

x 2 tan 2 0— y 2 K o M o'^- 


(71) 


This load distribution is shown in figure 18 for a triangular and a 
swept— back wing. It is seen that the loading at sonic speed bears 
a close resemblance to those found at higher Mach numbers. Two 
similarities of note are the discontinuity in the pressure gradient 
at the Mach wave originating from the trailing edge of the root 
chord nnri the satisfying of the Kutta condition only where the 
trailing edge is subsonic. The lift and induced drag coefficients 
are given, respectively, by 
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\ 


Cl = * a /, _ V 
a tan 0 2 tan 0 ^ s 0 2 . 



Figure 18. Pressure distributions for triangular and swept- 
back wings at M Q = 1. 

and 

C Dj = A 
a 2 tan 0 tan 0 

where 


k 2 tt 


4 (s 0 /c 0 tan 0) J 


(73) 



These coefficients are plotted as a function of A/tan 0 in figure 19. 
It is shown that the values of Cp^/a 2 tan 0 and Cj^/a tan 0 for 

finite aspect ratio swept-back wings are always less than the corre- 
sponding values for the triangular wing (A/tan 0=4). 
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When the span of the swept— hack wing becomes very large , the 

slope of the trailing edge 
approaches asymptotically the 
slope of the leading edge. It 
follows that for infinitely 
large aspect ratio the limiting 
value of the load distribution 
on the outboard sections should 
approach the value given by 
simple sweep theory for an 
infinitely long swept— back 
lifting surface with constant 
chord. This result is, in 
fact, a consequence of equa— 
tion ( 71 ). 
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Figure 19*— Variation of lift 

and drag with aspect ratio for 
a swept— back wing at M n = 1. 
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APPENDIX 


LIST OF IMPOKTANT SYMBOLS 


a i 

y coordinate of trailing edge 

a o 

free— stream speed of sound 

a 

local speed of sound 

a* 

critical speed of sound 


aspect ratio 


( span)" 


(wing area) _ 


NACA TN No. 182 ^ 


C(u) 

c o 

C L 


C Dl 


D 

D i 

erf (x) 


F n 

F Un> k n) 

K n 

V' V 

kn 


k* 


n 


hi 


Fresnel's cosine integral 
wing root chord 
lift coefficient 



u 


cos — x 2 dx 


) 


lift 


_q(wing area) 


dC L 

da 

indicial lift coefficient 

T i 

induced drag coefficient , — * r 

Lq(wing area) J 

drag 

induced drag 


x _^ 2 

error f unct ion of x ( — / e dX, 



elliptic integrals with moduli k' n 
modulus of elliptic functions 


l 


length of body 
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M o 

/V 0 \ 

free— stream Mach number 1 — J 

\ a o/ 

Ni(cut,v) 

C//oyt cQS _ sin x 2 s(o)t— x 2 ) dx 

«/v *- 

N 2 (<nt,v) 

r tjmt cos x 2 S(cct— x 2 ) + sin x 2 C(o)t— x 2 ) dx 

V7 L 

A£ 

q 

loading coefficient (pressure on lower surface 
minus pressure on upper surface divided by free— 
stream dynamic pressure) 

q 

free— stream dynamic pressure ^ 2 P o ^o 2 

r,e 

polar coordinates in yz plane (y = r cos 9, z = r sin0) 

s 

Fresnel's sine integral ^ j ^ sin ^ x 2 dx 

s 

operational equivalent of t 

s o 

wing semi span 

t * 

time 

• 

t 

a 0 t* 

t o 

maximum distance measured parallel to y axis from 
x axis to trailing edge (fig- 15) 

u,v,v 

perturbation velocity components in x,y, z directions, 
respectively 

V 

local velocity 

Vo 

free— stream velocity 


free— stream velocity 
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v r 

radial component of perturbation velocity 

x,y,z 

Cartesian coordinates 

a 

angle of attack in radians 

3 

y|l-M 0 2 ! 

7 

ratio of specific heats, for air 7 = 1.4 

A<p 0 ,Au 0 ,Aw 0 

discontinuity in component in z = 0 plane 

0 

semivertex angle of swept hack wing 

V 

CDCo 

2 

Po 

free— stream density 

<5 

total velocity potential 


perturbation velocity potential 

<D» 

impressed frequency (reference to true time) 

(0 

a o 
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